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The magnetization process of the isotropic Heisenberg antiferromagnet on the kagome lattice 
is studied. Data obtained from the numerical-diagonalization method are reexamined from the 
viewpoint of the derivative of the magnetization with respect to the magnetic field. We find that 
the behavior of the derivative at approximately one-third of the height of the magnetization 
saturation is markedly different from that for the cases of typical magnetization plateaux. The 
magnetization process of the kagome- lattice antiferromagnet reveals a new phenomenon, which 
we call the "magnetization ramp." 

KEYWORDS: kagome lattice, antiferromagnetic Heisenberg spin model, magnetization process, numerical- 
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Frustration has attracted considerable attention as an 
origin of various exotic phenomena in condensed-matter 
physics. In particular, one of the typical frustrated sys- 
tems in the field of magnetism is that due to the lattice 
structure based on the existence of neighboring bonds 
forming local triangles and tetragons. In such lattice 
structures, there is the case of the kagome lattice. 1 ) The 
5=1/2 kagome-lattice antiferromagnet has been stud- 
ied; although it is believed that no long-range order is 
realized even at zero temperature owing to the strong 
frustration and large quantum fluctuation, no conclusive 
evidence has been obtained so far in spite of extensive 
studies. Under these circumstances, the kagome-lattice 
antiferromagnet has recently become a hot topic again 
because of the discovery of some new materials. 2 6 ) 

From the theoretical point of view, however, it is 
well known that calculations of two-dimensional kagome- 
lattice systems are difficult even using a computational 
method. As reliable numerical methods, the quantum 
Monte Carlo (QMC) method, the density matrix renor- 
malization group (DMRG) method, and the numerical- 
diagonalization method are well known. Although the 
DMRG method can treat systems with large sizes, the 
dimensionality of the systems is limited to being less 
than two. QMC simulations can treat systems in higher 
dimensions, but the so-called negative sign problem pre- 
vents us from obtaining reliable results in the cases of 
frustrated systems. Nevertheless, double-peak behavior 
was clarified to appear in the specific heat as a result 
of the great effort devoted toward carrying out effective 
sampling in the QMC simulations. T ) Only the numerical- 
diagonalization method does not suffer from the limita- 
tion of dimensionality nor the negative sign problem. The 
disadvantage of the numerical-diagonalization method 
is the limitation that available system sizes are small. 
Thus, it is important to obtain precise results for sys- 
tems that are as large as possible in the available cal- 
culations by the numerical-diagonalization method and 
to examine the finite-size effect very carefully, which con- 
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tributes greatly the deep understanding of frustrated sys- 
tems. With this background, numerical-diagonalization 
studies of kagome-lattice systems 8-11 ) have been carried 
out. 

The magnetization process of the kagome-lattice 
Heisenberg antiferromagnet was examined by the 
numerical-diagonalization method. 12 15 ) In ref. 12, the 
full process for system sizes of up to N = 30 and the 
higher-field part of the process for N = 33 were reported. 
Hida pointed out that a magnetization plateau appears 
at one-third of the height of the saturation. In ref. 15, the 
magnetization process of the N = 36 system was given; 
the authors of ref. 15 claimed that the existence of the 
plateau is established from their result. 

In this letter, we reexamine the magnetization process 
of the kagome-lattice Heisenberg antiferromagnet from 
another viewpoint in analyzing finite-size calculations. 
The point is to observe the field-derivative of the magne- 
tization, namely, the differential magnetic susceptibility 
in magnetic fields. The observation will help us capture 
the behavior of the magnetization process with high ac- 
curacy. The purpose of this letter is to clarify, from the 
observation, that the magnetization process exhibits be- 
havior different from those reported so far, for example, 
the magnetization plateau and magnetization cusp. We 
will call the new behavior of the magnetization process a 
"magnetization ramp." For a system with size iV, we can 
obtain the lowest energy E(N, M) in each subspace for a 
given value of 5^ ot denoted by M by the numerical diag- 
onalization of the Lanczos algorithm and/or the house- 
holder algorithm, where Sl ot represents the z-component 
of the total spin. We can evaluate the derivative by the 
expression 



X 



E(N, M + 1) - 2E(N, M) + E(N, M+l) 



(1) 



1/M sat 

where M sat denotes the saturation of the magnetization, 
namely, M sat = NS for the spin-^ system. Note that x 
cannot be defined when one of the three E in eq. (1) does 
not become the ground-state energy of the system in any 
magnetic field. We obtain % in eq. (1) as a function of 
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Fig. 3. Shapes of the finite-size clusters in the kagome lattice. 
Cluster (C) is the same as that for N = 33 in ref. 12. Cluster 
(D) is the new cluster for N = 36, while cluster (E) is that for 
N = 36 in ref. 15. 



Fig. 1. Cluster shapes of the (A) one- and (B) two-dimensional 
systems of interacting S = 1 dimers. Our calculations have been 
carried out for N = 20 in both cases. In (B), the tilted square 
of y/20 x v^20 is shown with dotted lines. The thick green and 
thin black bonds denote the intra- and interdimer interactions, 
respectively. 
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Fig. 2. Field-derivative of the magnetization % of the one- and 
two-dimensional systems of interacting S = 1 dimers, denoted 
by crosses and circles, respectively. Insets (a) and (b) show the 
magnetization process of the one-dimensional case (A) in Fig. 1 
for J2/J1 = 0.15 and that of the two-dimensional case (B) for 
J2/J1 = 0.05, respectively. The main panel shows x as a function 
of the magnetization divided by the saturation with correspond- 
ing colors and symbols. 



M/M sat . 

Before investigating the kagome-lattice antiferromag- 
net, let us observe numerical-diagonalization finite-size 
data of x m one - and two-dimensional systems of inter- 
acting 5=1 dimers. This model reveals a typical mag- 
netization plateau at half the height of the saturation. 
The Hamiltonian of the interacting 5 = 1 dimers is given 

by 



H 



(2) 



where Si denotes the 5=1 spin operator. Here, the 
first term describes the intradimer interactions (denoted 
by the thick green bonds in Fig. 1), the second term de- 
scribes the interdimer interactions (denoted by the thin 
black bonds in Fig. 1), and the third term is the Zee- 
man term. Note that J\ is the energy unit; thus, we set 
Ji = 1. We examine the cases for the clusters depicted 
in Fig. 1. The results, which are depicted in Fig. 2, will 
help us accurately capture the behavior of x near the 
magnetization plateau. The system size is N = 20 com- 



monly. In the two-dimensional case, the tilted square 
shown in Fig. 1 is treated. The magnetization processes 
of the one-dimensional case (A) for J2/J1 = 0.15 and 
the two-dimensional case (B) for J2/J1 = 0.05 are pre- 
sented in insets (a) and (b) in Fig. 2, respectively. We 
choose these parameters so that the sum of the ampli- 
tudes of interdimer interactions operating a dimer of the 
thick green bond is the same in Figs. 1(A) and 1(B). We 
can clearly observe a plateau at M/M sat = 1/2 in both 
cases. Note that the experimentally observed magnetiza- 
tion processes were reported for one-dimensional systems 
of Ni compounds 16 ' 17 ) and for the two-dimensional sys- 
tem of the organic biradical F 2 PNNN0. 18 ) The one- and 
two-dimensional magnetization processes appear to be 
very similar; it is difficult to find a clear difference be- 
tween them. Let us next discuss the behavior of the 
results for which are presented in the main panel. Near 
half the height of the saturation, we can observe that x 
diverges in the one-dimensional system, while it remains 
finite in the two-dimensional system. These properties 
originate from the nature of the density of states de- 
termined from the parabolic dispersion. It is noticeable 
that the behavior does not change for each dimension, 
irrespective of the difference in the lower- or higher-field 
sides, and that the value of x just at M of the plateau 
is discontinuous with the values of x around it. One 
can find that x clearly exhibits the behavior of a typ- 
ical magnetization plateau from the above observation. 
Note that these characteristics of x appear not only at 
M/M sat = 1/2 but also at M = 0. The behavior at 
M = is already well known in various cases; 19 22 ) the 
behavior of x m these cases clearly satisfies the charac- 
teristics mentioned above. 

We now examine the magnetization process of the 
kagome-lattice antiferromagnet. Its Hamiltonian is given 

by 



n 



(3) 



where the sum of the first term runs over the nearest 
neighbors on the kagome lattice. Note that Si in eq. (3) 
denotes the 5 = 1/2 spin operator and that J is the en- 
ergy unit; thus, we set J = 1. The shape of the finite-size 
cluster of the kagome lattice under the periodic bound- 
ary condition is not necessarily unique even when the 
number of spin sites N is given. In fact, for TV = 36, 
there is another shape, denoted by the red cluster (D) in 
Fig. 3, which is different from that treated by Honecker 
et al. in ref. 15; the cluster of Honecker et al. is shown in 
Fig. 3(E). In this letter, we present the result for the red 
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Fig. 4. Results for the kagome- lattice antiferromagnets. (a) Mag- 
netization processes shown for the N = 36 cluster (D) and the 
N = 33 cluster (C) with circles and squares, respectively, (b) 
Field-derivative of the magnetization x as a function of the mag- 
netization divided by the saturation. Circles and squares in (b) 
denote the results of x corresponding to the cases in (a). Crosses 
correspond to the N = 36 cluster (E). Note that part of x f° r 
the TV = 36 cluster (E) is missing because x cannot be defined 
when the lowest-energy state in the subspace does not become 
the ground state of the system in the magnetic field. Diamonds 
represent the case of the N = 30 cluster investigated in ref. 12. 

cluster (D) for N = 36, which helps us know the finite- 
size effects from the difference between the two cases (D) 
and (E) for N = 36. 

The result of the magnetization process of the red clus- 
ter (D) for N = 36 sites is depicted in Fig. 4(a). We also 
give the full result of the blue cluster (C) for N = 33, 
which is the same cluster as that treated by Hida in 
ref. 12. In Fig. 4(b), we present the results of x f° r the 
finite-size systems of N = 30, 33, and 36. We focus our 
attention on the behavior of \ around M/M sat ~ 1/3. 
One can observe that around M/M sat ~ 1/3, shown in 
Fig. 4(b), x exhibits divergent behavior on the smaller-M 
side, while it is very small on the larger-M side; namely, 
the behavior on the smaller-M side is not the same as 
that on the larger-M side. This fact is markedly differ- 
ent from the behavior of x f° r the interacting S = 1 
dimer systems around M/M sat ~ 1/2 in Fig. 2. In par- 
ticular, the divergent behavior on the smaller-M side 
appears to be similar to that in the one-dimensional 
case even though the lattice is two-dimensional with- 
out any isotropy. This behavior may suggest that a one- 
dimensional spin structure forms in the two-dimensional 
lattice structure. A similar one-dimensional feature in the 
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Fig. 5. Schematic shape of the magnetization ramp. 

kagome lattice was reported in refs. 23 and 24. The rela- 
tionship between the present result and that in refs. 23 
and 24 is an open problem which should be tackled in 
the near future. Since the divergent behavior appears on 
the smaller-M side, the external field corresponding to 
M/M sat ~ 1/3 is the critical field. On the other hand, 
on the larger-M side, there is no divergent behavior of x- 
It is noticeable that x from the larger-M side is continu- 
ous with x j us t at M/M sat = 1/3. It is unclear whether 
or not x at M/M sat = 1/3 vanishes as N — >> oo at the 
present time. 

From these observations, the behavior of the magneti- 
zation around M/M sat = 1/3 is anomalous; it is reason- 
able to consider that the behavior is a new phenomenon 
in magnetization. The schematic behavior of M around 
M/M sat = 1/3 based on the discussion of the charac- 
teristics of x i s shown in Fig. 5. The shape is similar to 
the ramp of a Nordic ski jump when seen from the hor- 
izontally sideways direction. Thus, we call the behavior 
of M a magnetization ramp. This behavior is markedly 
different from the magnetization process of a classical 
system on the kagome lattice. 13 ) Although the magneti- 
zation ramp may appear to be a magnetization cusp, the 
ramp is definitely different from a cusp because a con- 
ventional cusp occurs at the level crossing between two 
states with finite gradients, for example, as reported in 
ref. 25. Thus, a magnetization ramp is an unusual and 
new phenomenon. 

Next, we discuss the characteristic behavior of other 
parts with respect to the relative height in our finite- 
size magnetization processes. We successfully observe a 
jump near the saturation, which was proven in ref. 26. 
Note that several eigenstates with various values of 5^ ot 
degenerate just at h/J = 3. In the N = 36 magneti- 
zation process reported in ref. 15, a jump appears at 
M/M sat = 11/18; at the field where the lowest-energy 
eigenvalue in the subspace of S'* ot = 10 meets that of 
gtot _ -L2 ? ^ e ener gy jg lower than the lowest-energy 
eigenvalue in the subspace of Sl ot = 11. In this meaning, 
this jump is different from the jump near the satura- 
tion discussed above. In the present result for TV = 36 in 
Fig. 4(a), however, the jump at M/M sat = 11/18 does 
not appear. It is thus unclear whether or not the jump at 
this relative height survives in the thermodynamic limit. 
Even in these situations, when one focuses on the behav- 
ior of x nea r M/M sat ~ 0.6 in Fig. 4(b), x 1S enhanced 
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irrespective of the system size except for the case of the 
green cluster (E) for N = 36, in which x cannot be ob- 
tained near M/M sat ~ 0.6; the enhancement suggests 
the existence of some anomaly around M/M sat ~ 0.6. It 
should be clarified in future studies whether this anomaly 
is a jump in the magnetization process or another. 

Finally, let us discuss the relationship between our 
observation of the magnetization process of the ideal 
S =1/2 kagome-lattice Heisenberg antiferromagnet and 
the magnetization measurement of the actual compounds 
volborthite and vesignieite, which have relatively small 
values of J. These materials include other factors that 
are not considered in the ideal situation assumed in this 
study. Vesignieite includes a small amount of impurity 
spins according to the results of recent studies. In vol- 
borthite, on the other hand, there is spatial anisotropy 
in the bonds, although very pure samples are available. It 
should be noted that the magnetization processes of both 
materials reveal flat behavior at about M/M sat ~ 0.4, 27 ) 
which is slightly larger than the (theoretical) numerical 
result of M/M sat = 1/3. It should be investigated in 
the near future whether or not the missing factors men- 
tioned above are the origin of this difference in height 
of M exhibiting the flat behavior and whether or not 
these factors affect the appearance of the magnetization 
ramp. In particular, in ref. 5 it was reported that in vol- 
borthite, another anomalous behavior in its magnetiza- 
tion process is observed at around M = (1/6)/ib/Cu and 
(1/45)/ib/Cu, referred to as "magnetization steps." 28 ) In 
the result for the red cluster (D) of TV = 36 in Fig. 4, 
small enhancements of x appear at M/M sat = 1/6 and 
1/18. Unfortunately, the resolution is not sufficient to 
judge whether or not the finite-size data surely capture 
the behavior of the magnetization steps; calculations for 
larger systems are greatly anticipated. 

In summary, we have examined the magnetization pro- 
cess of the S = 1/2 kagome-lattice Heisenberg antiferro- 
magnet from the viewpoint of its gradient x based on 
numerical-diagonalization finite-size data. The observa- 
tion of x is very useful for extracting characteristics in the 
magnetization process that are difficult to detect when 
only the magnetization process is observed. Our observa- 
tion leads to a new phenomenon, a magnetization ramp. 
To confirm the existence of the phenomenon, two ap- 
proaches will be necessary in future studies. One of them 
is the construction of an effective theory. The other is a 
numerical examination of systems with larger sizes. We 
plan to carry out calculations for N = 39; the results 
will be published elsewhere. In such works, something 
unknown in properties of kagome-lattice systems may be 
clarified further. 
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